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Abstract: The invariant spectra of 7r _ mesons produced in (ir + /K + )p interactions at 250 GeV/c are 
analysed in the framework of the hydrodynamical model of three-dimensionally expanding cylindrically 
symmetric finite systems. A satisfactory description of experimental data is achieved. The data favour 
the pattern according to which the hadron matter undergoes predominantly longitudinal expansion 
and non-relativistic transverse expansion with mean transverse velocity (u t ) = 0.20 ± 0.07, and is 
characterized by a large temperature inhomogeneity in the transverse direction: the extracted freeze- 
out temperature at the center of the tube and at the transverse rms radius are 140 ±3 MeV and 82 ±7 
MeV, respectively. The width of the (longitudinal) space-time rapidity distribution of the pion source 
is found to be A77 = 1.36 ±0.02. Combining this estimate with results of the Bose-Einstein correlation 
analysis in the same experiment, one extracts a mean freeze-out time of the source of (rf) = 1.4 ± 0.1 
fm/c and its transverse geometrical rms radius, Rc^ms) = 1.2 ± 0.2 fm. 
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1 Introduction 

Recent investigations of interference correlations of identical mesons (Bose-Einstein correlations, BEC) 
give evidence for a collective nature of the evolution of hadronic matter created in high-energy heavy 
ion collisions [1-4]. The hadronic matter flow (developed predominantly in the longitudinal direction) 
gives rise to dynamical correlations between space-time and phase-space coordinates of produced 
particles. As a consequence, the 'interferometric' radius, determining the shape of the two-particle 
correlation function at small relative momenta, turns out to be dependent on the kinematical variables 
of the particle pair. It does not measure the genuine geometrical size of the meson-emitting source, but 
the effective size of the source segment radiating mesons with sufficiently small relative momentum. 
The effective size reflects both the geometrical size and the 'homogeneity' length of the source. The 
latter is determined by macroscopic properties as temperature and expanding velocity profiles of the 
hadronic matter at freeze-out or hadronization time. 

These expectations are based on hydrodynamical models (see [5-7] and references therein) and are 
qualitatively confirmed by experimental data concerning interferometric radii and their dependence 
on the di-meson rapidity and transverse mass. It is interesting, that this observation not only holds 
for heavy-ion collisions [1-4], but also for meson-proton collisions [8,9] for which the formation of 
hadronic matter with properties of macroscopic (hydrodynamical) systems is less evident. 

One should stress, however, that correlation measurements do not contain the complete informa- 
tion on the geometrical and dynamical parameters characterizing the evolution of the hadronic matter. 
More comprehensive information can be provided by a combined analysis of data on two-particle cor- 
relations and single-particle inclusive spectra [5-7] . 

This work is, therefore, devoted to an analysis, in the framework of the hydrodynamical model 
for threc-dimcnsionally expanding systems [7], of the invariant spectra of tt~ mesons in the central 
rapidity region of (7r + /K + )p interactions at 250 GeV/c and its combination with earlier results [8,9] 
on two-particle correlations. The measurements are performed with the help of the European Hybrid 
Spectrometer (experiment NA22) at the CERN SPS. The first NA22 data on the pion inclusive 
spectra [10] (with about half the statistics used in the present paper) show their similarity in tt + and 
K + induced reactions. The present analysis, therefore, is based on the combined 7r + p and K + p data. 
Earlier results on BEC in the same data can be found in [8,9,11-14]. 

The hydrodynamical model parametrizations are discussed in Sect. 2. Experimental data and the 
results of their analysis are presented and combined with the results from two-particle correlations in 
Sect. 3. Conclusions arc summarized in Sect. 4. 



2 The parametrizations of invariant pion spectra 

The experimental data are analyzed in the framework of the hydrodynamical model for three-dimensio- 
nally expanding cylindrically symmetric systems [7] . There, the invariant spectrum of pions in rapidity 
y and transverse mass mt is approximated by 
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The width Ay of the rapidity distribution given by (2) is determined by the width Ar] of the 
longitudinal space-time rapidity rj distribution of the pion emitters and by the thermal smearing 
width ^/T /m t , where To is the freeze-out temperature (at the mean freeze-out time Tf) at the axis 
of the hydrodynamical tube, To = Tf(rx = 0). For the case of a slowly expanding system one expects 
Ar) <C To/mt, while for the case of a relativistic longitudinal expansion the geometrical extension Ar) 
can be much larger than the thermal smearing (provided m t > To). 

Except for the inhomogeneity caused by the longitudinal expansion, (1) also considers the inhomo- 
geneity related to the transverse expansion (with the mean radial component (itt) of hydrodynamical 
four-velocity) and to the transverse temperature inhomogeneity, characterized by the quantity 

(^)-^-i. (•) 

\ J- / J-rms 

where T rms = Tf(rx = rx(rms)) is the freeze-out temperature at the transverse rms radius rx(rms) 
and at time Tf. 

The exponential parameter a in (1) is related [7] to the number k of dimensions in which the 
expanding system is inhomogeneous. For the special case of the one-dimensional inhomogeneity (k = 
1) caused by the longitudinal expansion, a = 1 — 0.5fc = 0.5 (provided Arj 2 » To/mt). The transverse 
inhomogeneity of the system leads to smaller values of a. The minimum value of a = — 1 is achieved 
at k — 4 for the special case of a three-dimensionally expanding system with temporal change of local 
temperature during the particle emission process. 

The parameter yo in (1) denotes the midrapidity in the interaction c.m.s. and can slightly differ 
from due to different species of colliding particles. The parameter C is an overall normalization 
coefficient. 

Note that (1) yields the single-particle spectra of the core (the central part of the interaction 
that supposedly undergoes collective expansion). However, long-lived resonances also contribute to 
the single-particle spectra through their decay-products. Their contribution can be determined in 
the core- halo picture [20] by the momentum dependence of the intercept parameter \{y,mt) of the 
two-particle Bosc-Einstein correlation function. We determined this quantity in [8] in two different 
mt windows for the reaction investigated in the present paper, and found that within three standard 
deviations the A(mt) parameter is approximately independent of mt- Hence this correction can be 
absorbed in the overall normalization. 

The two-dimensional distribution (1) can be simplified for one-dimensional slices [7,15]: 
1. At fixed mt, the rapidity distribution reduces to the approximate parametrization 
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where C m is an mt-dependent normalization coefficient and yo is defined above. The width parameter 
Ay 2 extracted for different m t -slices is predicted to depend linearly on l/m t , with slope T and 
intercept Ar] 2 (cf. (2)). 

Note, that for static fireballs or spherically expanding shells (6) and (2) are satisfied with Ar] = 
[15]. Hence the experimental determination of the l/mt dependence of the Ay parameter can 
be utilized to distinguish between longitudinally expanding finite systems versus static fireballs or 
spherically expanding shells. 

2. At fixed y, the m 2 -distribution reduces to the approximate parametrization 

f(y,m t )=C y m?e X p(-^j (7) 



where C v is a y-dependent normalization coefficient and a is defined as above. 
The y-dependent "effective temperature" T c s(y) can be approximated as 



where T* is the maximum of T e ff(y) achieved at y — yo, and 

= T T, 
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with T and An 2 as defined above. 

The approximations (6) and (7) explicitly predict a specific narrowing of the rapidity and trans- 
verse mass spectra with increasing m t and y, respectively (cf. (2) and (8)). The character of these 
variations is expected [15] to be different for the various scenarios of hadron matter evolution. 

3 The results 

3.1 The data 

The experimental setup and the data-handling procedure are described in [10,16]. The results of this 
paper are based on the analysis of 104145 events of (7r + /K+)p interactions containing at least one 
negative track with momentum resolution better than 4% (depending on the momentum, the average 
momentum resolution varies from 1% to to 2.5%). All negative particles are assumed to have the 
pion mass. The contamination from other particles is estimated to be (7 ± 3)% [10]. Candidates for 
single diffraction dissociation with charge multiplicity n < 6 are excluded. For each event, a weight is 
introduced in order to normalize to the non-single diffractive topological cross sections [16]. 

The analysis is restricted to centrally produced 7r~-mesons with cms rapidity \y\ < 1.5. The total 
number of pions used in the analysis is about 397 k. 

3.2 The rapidity distribution 

The rapidity distributions fitted by (6) are given in Fig. 1 for 23 mt slices and mt values ranging 
from 0.16 to 0.63 GeV. The fit quality is satisfactory for all slices (0.7 < X 2 /NDF < 1.5). The fitted 
Ay 2 values plotted in Fig. 2 demonstrate the widening of the y-distribution with increasing l/m t . 
According to (2), the linear dependence of Ay 2 on l/mt is described by intercept Ar] 2 = 1.91 ± 0.12 
and slope To — 159 ± 38MeV. Thus, the width of the y-distribution is dominated by the spatial 
(longitudinal) distribution of pion emitters (inherent to longitudinally expanding systems) and not by 
the thermal properties of the hadron matter, as would be expected for static or radially expanding 
sources. 

3.3 The m t distribution 

The m t -distributions fitted by (7), are given in Fig. 3 for 11 y-slices for y-values between -1.5 to 
+1.5. For all slices x 2 /7VZ)F < 1.4. The fitted values of T e g(y) (Fig. 4) tend to decrease with 
increasing \y\ and approximately follow parametrization (8) with T» = 160 ±1 MeV, a = 0.083 ±0.007 
and yo = —0.065 ± 0.039. Note, however, that the fit of T e a(y) by (8) has low confidence level: 
X 2 /NDF = 29/11. This large value is due to an asymmetry in the T e g distribution with respect to 
y = 0. Except for the last point, T c g is higher in the meson than in the proton hemisphere. 

Note, furthermore, that the small value of the parameter a is the result of the comparatively large 
longitudinal extension An of the pion source (cf. (9)). Note also the consistency between the quoted 
value of a and the value a w 0.07 calculated from (9) at values of the parameters An 2 , To and T* 
extracted from two different fits of the one-dimensional spectra (6) and (7). 

The fitted values of the exponential parameter a in (7) are near zero (varying between — 0.02±0.02 
and 0.14 ± 0.02 for different y-slices). These small values of |a| correspond to a two-dimensional 
inhomogencity of the expanding system (a = 1 — 0.5fc). One concludes, therefore, that apart from 
a longitudinal inhomogeneity caused by the relativistic longitudinal flow, the hadron matter can 
also possess a transverse inhomogeneity (caused by transverse expansion or a transverse temperature 
gradient) or undergoes a temporal change of local temperature during the particle emission process. 



(9) 



3.4 The transverse direction 

Further information on hadron-matter evolution in the transverse direction can be extracted from a 
more detailed analysis using the two-dimensional parametrization (1) with parameters (u t ) and (^-) 
characterizing the strength of the transverse expansion and temperature inhomogeneity. The results 
of the fit are presented in Table 1. 

The fitted value of the exponential parameter of a = 0.26 ± 0.02 corresponds to an effective 
number k e s of dimensions in which the expanding system is inhomogeneous, k c s ~ 1.5, definitely 
in excess of k = 1 corresponding to a one-dimensional (longitudinal) inhomogeneity. This can be at 
least partly attributed to a transverse inhomogeneity of the system. The moderate value of the mean 
transverse four-velocity (u t ) = 0.20 ± 0.07 (Table 1) indicates that the transverse inhomogeneity is 
only to small extent caused by a (non-relativistic) transverse expansion. It is mainly stipulated by a 
rather large temperature inhomogeneity characterized by the fit result of ( A j£) = 0.71 ± 0.14. Using 
(5), one infers that the freeze-out temperature decreases from T = 140 ± 3 MeV at the central axis 
of the hydrodynamical tube to T rms = 82 ± 7 MeV at a radial distance equal to the transverse rms 
radius of the tube. 



3.5 Combination with two-particle correlations 

As already mentioned in the introduction, more comprehensive information on geometrical and dy- 
namical properties of the hadron matter evolution are expected from a combined consideration of 
two-particle correlations and single-particle inclusive spectra [5-7] . 

The two-particle correlation function K^ijPx,^) a ^ small momentum difference q = p x p 2 is 
often approximated by a Gaussian function [17]: 

K 2 ( Pl , p 2 ) ~ eM-Rhl - #o?out - ^s<4icL (io) 

where the three orthogonal components , q out , <7 s idc of the vector q arc oriented, respectively, along 
the collision axis, along and perpendicular to the pair transverse momentum; Rl,R and R s are, 
respectively, the longitudinal, 'out' and 'side' effective dimensions of the source segment radiating the 
BE correlated pion pairs. Due to the non-static nature of the source, these effective sizes vary with 
the average transverse mass M t — \ {m tl + m t 2) and the average rapidity Y = \{yi + V2) of the pion 
pair. In the 'longitudinal c.m.s.' (LCMS) [18], where Y = 0, the effective radii can be approximately 
expressed as [7,15,19]: 

Rl=T?Arit (11) 
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where parameters An 2 , T , (u t ) and (-^f) are defined and estimated above from the invariant spectra; 
Rq is related to the transverse geometrical rms radius of the source as i?G( rm s) ; Tf is the 

mean freeze-out (hadronization) time; At* is related to the duration time Ar of pion emission and to 
the temporal inhomogeneity of the local temperature; if the latter has a small strength (as one can 
deduce from the restricted inhomogeneity dimension estimated above: k c g « 1.5), an approximate 
relation At > At, holds; the variable (3 t is the transverse velocity of the pion pair. 

Relations (11)— (15) , combining the results of the single-particle invariant spectrum and the 
two-particle correlation function, allow one to extract additional information on the parameters 
(ti , At*, Rq) characterizing the space-time evolution of hadron matter. 

The interferometric radii Rl,R 07 R s were extracted recently [8,9] in the same experiment from 
the correlations of negative pions. To match with parametrization (10), the radii quoted in [8,9] are 
divided here by a factor 



The effective longitudinal radius i?L, extracted for two different mass ranges, M t = 0.26 ± 0.05 
and 0.45 ± 0.09 GeV/c 2 arc found to be R L = 0.93 ± 0.04 and 0.70 ± 0.09 fm, respectively. This 
dependence on M t matches well the predicted one. Using (11) and (14) with To = 140 ± 3 MeV and 
An 2 = 1.85 ±0.04 (Table 1), one finds that the values of Tf extracted for the two different Mt-regions 
are similar to each other: Tf = 1.44 ± 0.12 and 1.36 ± 0.23 fm/c. The averaged value of the mean 
freeze-out time is Tf = 1.4 ± 0.1 fm/c. 

Since we find that Arj is significantly bigger than 0, static fireballs or spherically expanding shells, 
that were found to be able to describe our two-particle correlation data in [8], fail to reproduce our 
single-particle spectra. 

The transverse-plane radii R and R s measured in [8,9] for the whole Mt range are: R = 0.91 ± 
0.08 fm and R s = 0.54 ± 0.07 fm. Substituting in (12) and (13), one obtains (at (i t = 0.484c [8]): 
At* = 1.3 ± 0.3 fm/c. The mean duration time of pion emission can be estimated as At > At, = 
1.3 ±0.3 fm/c. A possible interpretation of At « Tf might be that the radiation process occurs during 
almost all the hydrodynamical evolution of the hadronic matter produced in meson-proton collisions. 

An estimation for the parameter Rq can be obtained from (13) and (15) using the quoted values 
of R s ,To, (u t ) and (-^?) at the mean value of (M t ) = 0.31 ± 0.04 GcV/c (averaged over the whole 
Mf-range): Rg — 0.88 ±0.13 fm. The geometrical rms transverse radius of the hydrodynamical tube, 
i?G( rms ) = V%Rg = 1-2 ± 0.2 fm, turns out to be larger than the proton rms transverse radius. 



3. 6 Estimation of systematic errors 

The analytic formulae used to estimate the hydrodynamical model parameters were evaluated under 
certain approximations. As shown in [21], an accuracy better than 10-20% is reached under the 
following conditions: 



T R G T /M t 



With our parameter values, the l.h.s. of (16) turns out to be less than 0.13 for the whole range of M t ; 
b) Additionaly, the validity of the formulae (l)-(9) for the invariant spectra requires: 

\y-yo\ <1 _ (17) 



1 + A77 2 m t /To 

With our parameter values, the l.h.s. of (17) turns out to be less than 0.52 for the full ranges of 
y and m t considered. 

c) Formulae (10)-(15) for the two-particle correlation function are derived under the additional, 
more stringent condition: 

lY - yo1 <1. (18) 



l + A77 2 (M t /T -l) 

At the maximum value of \Y — y a \, (18) is fulfilled, except for a fraction of 8% of pion pairs 
with M t < 0.18 GeV/c 2 . In the full range of \Y — yo\ < 1.5 considered, this is reduced to less than 
4% of pion pairs. Note further that the upper limit of 10-20% relative errors is reached only at the 
upper limit of inequalities (16-18) and that the squared relative errors increase quadratically with the 
increase of the left-hand-sides of (16-18). 

Hence we estimate an upper limit of order 10-20% x (0.13/0.6) x (0.52/1.0) x 1.0 =1-2% for the 
precision of the analytic formulae used in the fits, except for 4% of the pion pairs where the relative 
systematic errors may reach 10-20%. This accuracy is well within the mean experimental (statistical) 
errors of the measured two-dimensional invariant spectra, analyzed in the present work (about 10%) 
and of the three-dimensional correlation function (about 15%), analyzed in [8,9]. 



3.7 The space-time distribution of ir emission 

On Figure 5 we show a reconstruction of the space-time distribution of pion emission points, expressed 
as a function of the cms time variable t and the cms longitudinal coordinate z. The momentum- 



integrated emission function along the z-axis, i.e., at r t = {r x ,r y ) = (0,0) is given by 

S(t, z) oc exp {-^-r- ) exp ^—^5- ) • (19) 

It relates the parameters fitted to our data with particle production in space-time. Note that the 
coordinates (t, z), can be expressed with the help of the longitudinal proper-time r and space-time 
rapidity r\ as (t cosh(?7), t sinh^)). 

We find a structure looking like a boomerang, i.e., particle production takes place close to the 
regions of z = t and z = —t, with gradually decreasing probability for ever larger values of space- 
time rapidity. Although the mean proper-time for particle production is Tf = 1.4 fm/c, and the 
dispersion of particle production in space-time rapidity is rather small, Ai] = 1.35 fm, we still see 
a characteristic long tail of particle emission on both sides of the light-cone, giving a total of 40 fm 
maximal longitudinal extension in z and a maximum of about 20 fm/c duration of particle production 
in the time variable t. 

4 Summary 

The invariant spectra of 7r~-mesons produced in (7r + /K + )p-interactions at 250 GeV/c are analyzed 
in the framework of the hydrodynamical model of a three-dimensionally expanding cylindrically sym- 
metric finite system [7]. 

The data favour a picture according to which the hadron matter undergoes extensive longitudinal 
expansion with a space-time rapidity width of An = 1.36 ± 0.02 and a non-relativistic transverse 
expansion with the mean transverse four-velocity (u t ) — 0.20 ± 0.07. 

The hadron matter possesses a large temperature inhomogeneity in the transverse direction: the 
freeze-out temperature varies from T = 140 ± 3 MeV at the central axis of the hydrodynamical tube 
to T rms = 82 ± 7 MeV at a radial distance equal to the transverse rms radius of the tube. 

The information from the single-particle invariant spectrum is combined with that from the two- 
particle correlation function. The transverse mass dependence of the width Ay of the single-particle 
y-distribution and that of the interfcrometric longitudinal radius i?L are found to be consistent. This 
allows to extract the mean freeze-out (hadronization) time of the strongly interacting matter: Tf = 
1.4 ± 0.1 fm/c. The duration time of pion emission is estimated to be At > 1.3 ± 0.3 fm/c, i.e. close 
to Tf, indicating that the emission process occurs during almost all the hydrodynamical evolution. 
As an estimate for the transverse geometrical rms radius of the hydrodynamical tube we obtain : 
i? G (rms) = 1.2 ±0.2 fm. 

The combined analysis of single-particle spectra and Bose-Einstein correlation functions increases 
the selective power of the analysis. Static (Kopylov-Podgoretskii) fireballs and spherically expanding 
shells, that were found to be able to reproduce our interferometric data in [8] fail to reproduce 
simultaneously the single-particle spectra and the BE correlations in the NA22 experiment. 

The systematic errors on the quoted parameters related to the limited accuracy of the approximate 
analytic formulae (used in the present study) are estimated to be of the order of or smaller than the 
quoted statistical errors. 

As far as we know, this is the first time that a full reconstruction of the space-time distribution 
of the particle emitting source is obtained in hadron-hadron reactions from a combined analysis of 
single-particle spectra and Bose-Einstein correlation measurements. 
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Table 1: Fit results according to parametrization (1) for \y\ < 1.5. 



Figure Captions 



Fig. 1 The rapidity distributions of centrally produced pions (\y\ < 1.5) for different m t -slices given. 
The curves are the fit results obtained according to parameterization (6). 

Fig. 2 The (l/m t )-dcpcndence of (Ay) 2 for inclusive ir~ meson rapidity distributions at \y\ < 1.5. The 
straight line is the fit result according to parametrization (2). 

Fig. 3 The mt-distribution of centrally produced pions for different y-slices, as indicated. The solid 
lines are the fit results obtained according to parameterization (7). 

Fig. 4 T e ff as a function of y fitted according to parametrization (8). 

Fig. 5 The reconstructed S(t, z) emission function in arbitrary units, as a function of time t and 
longitudinal coordinate z. The best fit parameters of Arj = 1.36, yo = 0.082, At = 1.3 fm/c 
and Tf = 1.4 fm/c are used to obtain this plot. 
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